private double IntFun(TVec x, TVec c, params object[] o) { double x = x[0]; return System.Math.Sin(x)*System.Math.Exp(-x*x); } private void Example() { Vector bpoints = new Vector(0); Vector weights = new Vector(0); MtxIntDiff.WeightsGauss(10,bpoints,weights); double area = MtxIntDif.QuadGauss(IntFun,-0.5*System.Math.PI,System.Math.PI,bpoints,weights,64); }